Metastable 7r-junction between an s±-wave and an s-wave superconductor 
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We examine a contact between a superconductor whose order parameter changes sign across 
the Brillioun zone, and an ordinary, uniform-sign superconductor. Within a Ginzburg-Landau type 
model, we find that if the the barrier between the two superconductors is not too high, the frustration 
of the Josephson coupling between different portions of the Fermi surface across the contact can 
lead to surprising consequences. These include time-reversal symmetry breaking at the interface 
and unusual energy-phase relations with multiple local minima. We propose this mechanism as a 
possible explanation for the half-integer flux quantum transitions in composite niobium-iron pnictide 
superconducting loops, which were discovered in a recent experiment pQ. 



Introduction- Understanding the structure of the order 
parameter of the iron-based pnictide superconductors [2] 
is the key to unveil their pairing mechanism. A conven- 
tional, phonon-mediated mechanism is usually associated 
with an order parameter of a uniform sign, while sign 
changes in the order parameter are typical of unconven- 
tional mechanisms, in which pairing is driven by purely 
repulsive Coulomb interactions. 

In the cuprates, the unambiguous identification of the 
d-wave symmetry of the order parameter came primar- 
ily from phase-sensitive experiments [3[ [4]. These experi- 
ments exploit the fact that due to the d-wave symmetry, 
specific geometries (such as a corner junction with an 
s-wave superconductor, or a tri-junction between three 
d-wave superconductors) are guaranteed to produce a tt 
phase shift in the phase of the superconducting order pa- 
rameter. In the iron arsenides, similar experiments [5j |6] 
have found no evidence for d-wave symmetry. 

From the theory side, it has been proposed that the 
pnictides have an extended s-wave ("s± wave") order 
parameter [3 [8], which can be chosen to be real but 
changes sign between the electron and hole pockets, and 
is invariant under the overall tetragonal symmetry (Ai p ). 
Such a pairing state has been found from solving both 
weakly [9UT2 and strongly [T3UT5] interacting models. 
Although there exist several experimental indications 
that this is indeed the correct pairing state in certain 
pnictide superconductors [T6UT9] . more direct experimen- 
tal evidence is highly desirable. Designing a Josephson 
interferometry device which could detect the s± state 
poses a significant challenge, since symmetry alone does 
not guarantee a tt phase shift in any geometry. Several 
ideas have been proposed to overcome this difficulty [20]- 
[24] . but none were realized to date. 

Progress has been made recently, in the work of Chen 
et al.[l . In this experiment, the flux through a com- 
posite Niobium(Nb)-NdFeAsOo.88Fo.i2(FeAs) supercon- 
ducting loop was measured. By using an external elec- 
tromagnetic pulse, is was shown that both integer and 



half-integer flux jumps can be induced in the loop, in 
units of the superconducting flux quantum <3>o = hc/2e. 
While providing a strong indication for a sign change in 
the order parameter, these results are surprising, because 
neither a tt-j unction nor a 0-j unction between the Nb and 
the FeAs superconductors would lead to the possibility 
of half-integer flux quantum jumps. 

Motivated by these experiments, we study a model of 
a junction between a sign-changing and a conventional 
(s-wave) superconductor. The model allows us to inter- 
polate between the tunneling (weak coupling) and the 
metallic contact (strong coupling) regime. We find that 
above a certain critical coupling strength, the frustration 
of the Josephson coupling across the barrier can lead to 
unusual energy-phase relations in the junction. If the 
phase stiffness of the superconductors is small enough, 
the energy-phase relation has two minima which break 
time-reversal symmetry; if the phase stiffness is large, 
the global minimum of the energy is at a phase differ- 
ence of Aip = 0, but an additional met a- stable minimum 
at Aip = 7r appears. The latter situation can explain 
the Chen et al. experiment, since it allows allows the 
junction to switch from Aip = to 7r, causing a half- flux 
quantum jump in the loop. The fact that the half-integer 
jumps appear only beyond a certain value of the critical 
current in the loop, as well as the relatively small prob- 
ability of half-integer jumpsfl , can both be understood 
within our model. 

The model- We consider a planar Josephson junction 
between an s-wave and a sign- changing s-wave, shown in 
Fig. [I] Since we are interested in the crossover from 
a tunneling barrier to a metallic contact, in which the 
order parameters of both superconductors are modified 
significantly at distances of the order of a few coher- 
ence lengths from the junction [25], the problem needs 
to be treated self-consistently. Rather than solving the 
full Bugoliubov-de Gennes equations, we use an effective 
Ginzburg-Landau type free energy functional which de- 
pends on the superconducting order parameter near the 



junction [26] . Although the Ginzburg-Landau descrip- 
tion is strictly valid only close to the critical temperature 
(T c ) of both superconductors, we expect it to capture the 
qualitative behavior of the system even at lower temper- 
atures. 

In order to capture the multi-band nature of the sys- 
tem, we introduce two superconducting order parame- 
ters, A^ where i = 1,2. Microscopically, these can be 
viewed as belonging to regions of different momenta par- 
allel to the junction: 
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FIG. 1: (Color online.) Junction between an s- wave and an s± 
wave superconductor. We define two order parameters, Ai=i^ 
(see Eq. [I]), which belong to different regions in momentum 
space. In the s± side, Ai,2 have an opposite sign. 



where ip k a (x) is the electron creation operator at po- 
sition x, momentum fcy parallel to the junction, and 
spin cr, Vfe^fe/ is the pairing interaction in the Cooper 
channel, x is the coordinate perpendicular to the junc- 
tion, A is the area of the junction, and the two momen- 
tum regions are defined by Qi — j^y fco > |^|||| and 

^2 = | fcy |fco < l^lll}' wnere is an arbitrary momen- 
tum chosen such that in the s± side, Ai > and A 2 < 0. 
There, Ai and A 2 can be thought of as the order param- 
eters on different bands. Note that such as a decomposi- 
tion is possible irrespective of the relative orientation of 
the two crystals, as long as in the sign-changing s-wave 
side, the region £l\ (Q 2 ) is dominated by the positive 
(negative) part of the order parameter. 

We describe the junction using the following phe- 
nomenological Ginzburg-Landau free energy: 



F[Ai,A 2 ] =F L + F R + F C 



(2) 



Here, 



F V [A U A 2 ] = J^dxl^ ^\ dxAi \ 2 -^\Ai 



~< |A,| 4 + I< (A*d x A z - A^A*) 2 



-w"(A*A 2 + c.c.) ), 



(3) 



where F v= l^r are the free energies of the left (Nb) and 
right (FeAs) sides, with J L = f® £ and J R = J* (21 is the 
system length), and 



F c [Ai, A 2 ] = [ T * I A * ( 0+ ) - A * (°~ 

i=l,2 



(4) 



describes the contact between the two superconduc- 
tors. In some cases, we will consider a "ring" 
geometry in which we add to F c another term, 

2" 
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= 1,2 



T;| A; (£)-Ai (-£)r 



k\ ,r% in Eq. [3] are the standard Ginzburg-Landau 
parameters of band i = 1,2 on the left /right (y = L, R 
respectively) of the i unction [27]. The ol\ terms repre- 
sent an additional energy cost of creating super currents. 
These terms turn out to be necessary to describe the 
transition from a single-minimum to a double-minimum 
junction (see below). Close to the critical temperature, 
the a term becomes negligible compared to the other 
terms in Eq. [3J since it is of higher order in the A^'s and 
their derivatives. At lower temperatures, however, it can 
become important. 

The parameter v v describes the inter-band coupling, 
and encodes the tendency towards s or s± pairing: pos- 
itive (negative) v L,R corresponds to an s- (s±-) wave 
superconductor, respectively. 

Ti (Eq. [4| represent the strengths of the couplings of 
the two order parameters across the barrier. These pa- 
rameters allow us to interpolate between the tunneling 
regime (Ti —> 0) to the metallic contact regime (Ti — » oo). 
Note that in the latter regime, Ai j2 become continuous at 
the junction. This is the boundary condition of a metal- 
lic contact between two superconductors, assuming for 
simplicity that the density of states at the Fermi energy 
times the pairing interaction is continuous at the con- 
tact [25] . Our results do not depend qualitatively on this 
assumption. 

Results- We minimize the free energy (Eq. |5| numer- 
ically. The minimization is done on a discrete lattice, 
with a small enough lattice spacing such that the results 
are independent of its size. In order to map the energy- 
phase relation of the s-s± junction, we use the follow- 
ing boundary conditions at x = ±1: argAi(— £) = 0, 
argAi(^) = A(p, where cp is varied between and n. 
The results (in particular, the qualitative behavior of the 
junction) do not depend on the choice of I. 

Depending on the contact strength and on various ma- 
terial parameters, we find three qualitatively different 
regimes: 1. A "single minimum" regime, which is re- 
alized for small T^, in which the free energy is minimal 
at Acp = 2n7T (where n is an integer); 2. A "time rever- 
sal breaking" (TRB) regime [28] at intermediate Ti and 
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FIG. 2: (Color online.) (a) Phase diagram of an s — s± junc- 
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tion, using = 1, k\ 
ol\ — a, £ = 3 and Ti —T^. The different states correspond to 
different qualitative behaviors of the energy-phase relation in 
the junction: a "single minimum" state in which the minima 
of the free energy occur at a phase difference of 2n7r, a "time- 
reversal breaking" (TRB) phase in which the minima occur 
at 2n7r± Acpo, and a "double minimum" phase in which there 
are minima at mr. The triangle, the circle and the square 
correspond to the three parameter sets which are used in Fig. 
[3] (b) Same as (a), for a case where the parameters are tuned 
such that Ai,2 are equivalent (by setting u\ — u% = 1). In 
this case, the single minimum phase does not occur. 



small a, where the free energy exhibits degenerate min- 
ima at Acp = 2nir ± cpo; 3. A "double minimum" regime, 
in which there are global minima at Acp = 2nir and local 
minima at Acp = (2n + 1)tt. This regime is realized for 
sufficiently large a and Ti. 

The phase diagram as a function of a and T i: showing 
the three phases described above, is shown in Fig. [2^i. 
The parameters used are listed in the figure caption. 

In order to understand these results, we consider an 
artificial situation in which the parameters of the free en- 
ergy are tuned such that the two order parameters Ai 5 2 
are exactly equivalent. In this case, the free energy of the 
junction F(Acp) (after minimization over Ai^) is invari- 
ant under a shift of Acp by tt. Together with time-reversal 
symmetry, F(Acp) = F(—Aip), this dictates that there 
are two generic situations: F(Acp) is minimal for either 
Acp = nir or Acp = (n + 1/2)tt. This can be understood 
as follows: due to the sign change of the order parameter 
in the FeAs side, the Josephson coupling is frustrated. 
In order to relieve this frustration, the system can either 
twist the relative phase of Ai^ close to the interface, or 
deform the amplitudes of Ai^ such that |Ai| ^ | A.2 1 - 
The first effect favors Acp = (n + 1/2) 7r, and the second 
favors Acp = utt. Which effect dominates depends on the 
energy cost of a non-uniform phase relative to that of a 
non-uniform amplitude of the order parameter. Increas- 
ing a/ft increases the phase stiffness, and hence drives a 
transition between the two regimes. The phase diagram 
for the case of equivalent A 12 is shown in Fig. [2jo. 

Upon making A\^ inequivalent, the period of F(Acp) 
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FIG. 3: (Color online.) Free energy as a function of flux in a 
composite s-s± loop (inset) of length 2£ = 20, for Ti,2 — 2.67 
and ol/k — 0,0.25,1.25. These correspond to the "time- 
reversal breaking" , "single minimum" and "double minimum" 
states, respectively. The self-inductance of the loop was cho- 



sen such that H/2L = 0.1. 



becomes 27r. However, the local minima of F(Acp) remain 
stable over a finite region in parameter space. Therefore, 
the two phases described above survive. Between them, 
a third phase with minima at Acp — 2nir appears. 

A composite s-s± superconducting loop (shown 
schematically in Fig. [3| is also described by three qual- 
itatively different regimes, corresponding to the three 
phases of a single junction. The free energy of the loop 
as a function of the flux through the loop is given by [29] 



F loop = 



$2 
2L' 



(5) 



where F(J&) is the free energy of Eq. [2] with d x 
d x — 2eA/h (A = &/2£ is the vector potential), supple- 
mented with an additional contact term (of the form of 
Eq. [4] with x = tj , and L is the self-inductance of the 
loop. Fig. [3] shows the free energy of the loop as a func- 
tion of <£, with parameters (a,T) which place the system 
in one of the three distinct regimes described above. For 
intermediate values of a/ ft, the free energy of the loop has 
minima at approximately = n^o where n is an integer 
("single minimum" regime). For smaller a/ft, each mini- 
mum splits into two degenerate minima at <£> = n^oiA^ 
where A<£ depends continuously on the various parame- 
ters of the junction (TRB regime). Finally, for a/ft larger 
than a critical value, minima appear both at $ = n^o 
and <fr = (n + 1/2) $0 ("double minimum" regime). 

Comparison to the Chen et al. experiment- The ex- 
periment of Chen et a/.[l can be interpreted in terms of 
our model as follows. As the critical current in the com- 
posite Nb-FeAs loop is increased, the Nb-FeAs junction 
undergoes a transition from a single minimum to a dou- 
ble minimum phase. Upon entry to the double minimum 
phase, both n<3>o and (71+ 1/2) $0 flux jumps are observed, 
and the probability for (n + l/2)$o jumps increases with 
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increasing critical current (which corresponds to increas- 
ing Ti >2 in Fig.§. 

Whenever the composite loop is excited electromag- 
netically, the flux through the loop acquires a random 
"kick" , and the system can jump from one local minimum 
to another. The probability of ending at a metastable lo- 
cal minimum, <3> = (n + l/2)<I>o, is smaller than the prob- 
ability of ending in a global minimum, = n$o- There- 
fore, if the system jumps from ni<£>o to (77,2 + 1/2) $0, the 
next jump is likely to be to n^^o (with integer m, ri2, 
77,3). This implies that the half flux quantum jumps are 
correlated, and tend to appear in pairs. Such a trend 
is clearly visible in the experimental results [TJ [30] . We 
conclude that the experimental results of Ref. [1 can be 
naturally explained by the appearance of a metastable 
local minimum in the energy-phase relation of the Nb- 
FeAs junction at Aip = (n + l/2)7r, corresponding to the 
"double minimum" phase in our model. 

Conclusions- We have presented a model for a Joseph- 
son junction between a simple s-wave and a sign- 
changing s± wave superconductor. Due to the sign 
change in the s± side, the Josephson coupling across the 
junction is partially frustrated. If the barrier between 
the two superconductors is low enough, corresponding to 
a metallic contact, the system can relieve this frustra- 
tion by either breaking time-reversal symmetry, or de- 
veloping additional local minima at a phase difference 
of (2n + 1)tt. The half-flux quantum jumps in the ex- 
periment of Chen et al. can be explained by the second 
scenario. (Note that a time-reversal breaking junction 
would correspond to fractional flux jumps which are nei- 
ther n$o nor (n + 1/2) $0, which were not observed.) 

It is important to note that this behavior is unique to 
a junction between an s-wave and a sign changing super- 
conductor. In a junction between two s-wave supercon- 
ductors, only the single minimum phase is realized, for 
any strength of the coupling T\^- To verify this, we con- 
sidered a case in which v L = v R > 0, i.e. both sides of 
the junction are s-wave superconductors, and found only 
a single- minimum phase. In this respect, the observation 
of half- flux quantum jumps is a strong indication of a 
sign-changing order parameter. Together with the lack 
of observation of spontaneous flux in a polycrystalline 
sample [5 , which essentially rules out a d-wave order pa- 
rameter, we conclude that the s± is the most likely can- 
didate for the order parameter of F-doped NdFeAsO. 

The presence local minima in the energy-phase relation 
of a contact between a conventional and an iron-based su- 
perconductor can be looked for in experiments. The AC 
Josephson effect should reveal pronounced high harmon- 
ics of the AC current as a function of the bias voltage. In 
a composite loop, of the same kind as the one studied by 
Chen et al, half- flux quantum entries should be observed 
as a function of an external field. 

On the theoretical side, it will be interesting to ex- 
amine a strongly-coupled s-s± junction in a microscopic 



model. Within our phenomenological model, we found 
that the fate of the junction in the strong coupling limit 
is determined by the ratio of the phase and amplitude 
stiffnesses. In a Bardeen-Cooper-Schrieffer superconduc- 
tor at low temperatures, these stiffnesses both scale in 
the same way - their magnitude is proportional to the 
Fermi energy. Therefore, the relative magnitude of the 
two stiffnesses is sensitive to microscopic details. It will 
be particularly interesting to look for parameter regimes 
in which the broken time-reversal junction can arise. 
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I 

SUPPLEMENTARY INFORMATION: PAIRING OF HALF-FLUX JUMP EVENTS 

In Fig. [I] we show scans of the flux vs. time from the Chen et. al. experiment pQ . The arrows indicate the position 
of fractional (n + 1/2) 3>o flux jumps. From the data, it is apparent that the (n + 1/2) 3>o jumps are paired, i.e., a 
half-integer flux jump tends to be followed by another half-integer jump. This is consistent with the picture of a local 
minimum in the energy vs. flux relation of the loop at (n + l/2)$o- After each jump, the flux is less likely to get 
trapped in a local minimum than in a global minimum at n$o- Therefore, a jump from a global to a local minimum 
(which is accompanied by a change in the flux by (n + 1/2) $0) is likely to be followed by a jump back to the global 
minimum, which also changes the flux by (n + l/2)$o- ft would be useful to perform a more careful statistical analysis 
of the data, which could quantify the degree of correlation of half flux events, and possibly shed more light on the 
nature of this process. 
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FIG. 4: Flux vs. time in a composite Nb-NdFeAsO0.8sF0.12 from Chen et. al.pQ. The abrupt jumps in the flux are caused 
by external electromagnetic pulses. In each jump, the flux changes by either n$o or (n + l/2)<3>o, where n is an integer and 
$o = hc/2e. The arrows indicate the position of the (n + l/2)$o jumps. 



